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Abstract. We investigate the problem of controlling the magnetic moment in a 
ferromagnetic nanowire submitted to an external magnetic field in the direction 
of the nanowire. The system is modeled with the one dimensional Landau- 
Lifschitz equation. In the absence of control, there exist particular solutions, 
which happen to be relevant for practical issues, called travelling walls. In this 
paper, we prove that it is possible to move from a given travelling wall profile 
to any other one, by acting on the external magnetic field. The control laws 
are simple and explicit, and the resulting trajectories are shown to be stable. 

1. Introduction and main result 

The most common model used to describe the behavior of ferromagnetic mate- 
rials, called micromagnetism, was introduced by W.-F. Brown in the 60's (see T). 
It is based on a thermodynamic approach, and the main idea is that equilibrium 
states of the magnetization minimize a given energy functional, consisting of several 
components. The main components, which permit an accurate description of the 
behavior of ferromagnetic materials, are the magnetostatic one (electromagnetism) , 
the exchange one (spin-spin interactions), the anisotropy one (crystal shape in- 
fluence) and the Zeeman one (external influences). This point of view permits to 
recover the standard dynamical approach of ferromagnetic phenomena, based on the 
so-called Landau- Lifschitz equation, which was introduced in the 30's in 16 . This 
equation contains a haniiltonian term stemming from the Larmor microscopic spin 
precession equation, and a purely dissipative term, perpendicular to the precession 
component and related to the Euler equation of the static energy functional. 

More precisely, ferromagnetic materials are characterized by a spontaneous mag- 
netization described by the magnetic moment u which is a unitary vector field link- 
ing the magnetic induction B with the magnetic field H by the relation B = H + u. 
The magnetic moment u is solution of the Landau-Lifschitz equation 

(1) — = -u/\H,-u/\{uAH^), 

where the effective field is given by He = Au -I- hd{u) + Ha, and the demagnetizing 
field hd{u) is solution of the magnetostatic equations 

div B = div (H + u) ^0 and curl H = 0, 
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where Ha is an applied magnetic field. More details on the ferromagnetism model 
are provided in ^2 El Existence results have been established for the 

Landau-Lifschitz equation in ^HHIEHSOIj numerical aspects have been investigated 
in ^3 El , and asymptotic properties have been proved in El El ■ 

In this article, we consider an asymptotic one dimensional model of ferromagnetic 
nanowire submitted to an applied field along the axis of the wire. Let (61,62,63) 
denote the canonical basis of IR^ . The ferromagnetic nanowire is assumed to have 
an infinite length, and is represented by the axis Mei. The demagnetizing energy is 
given by hd{u) = —^262 — ^363 where u — (mi,U2,M3) (see ^21 where this formula 
has been derived using a BKW method, by considering a wire of nonzero diameter, 
and taking the limit when the diameter tends to zero). In addition, we assume that 
an external magnetic field S{t)ei is applied along the wire axis. The real-valued 
function S{-) is our control. 

The Landau-Lifschitz equation writes 

(2) — ^ -uAhs{u) ~uA{uAhs{u)), 

where hs{u) = |^ - U262 - ^363 -I- Sei. 

The magnetic field S{-)ei is generated by a coil winding up around the nanowire. 
Note that it would be extremely costly to generate a magnetic field in other direc- 
tions, all along the nanowire. By the way, such other controls do not lead actually 
to relevant solutions. 

In terms of control system, setting h(u) — u^x — ^262 — U363, Equation jSJ yields 

(3) Uf = -M A h{u) - u A (u A h{u)) - S{u A ei + u A {u A 61)). 

For 6 = 0, physical experiments demonstrate the existence of a particular sta- 
tionary solution, splitting the nanowire into two parts. The magnetic moment is 
almost equal to 61 in one of them, and to —61 in the other. This fundamental 
stationary solution, called a wall, is analytically given by 



Here, and throughout the paper, the notations ch , sh , and th , respectely stand for 
the hyperbolic cosine, sine, and tangent functions. 

Stability properties of the solution AIq for System Q with ^ = have been 
established in 0. 

When applying a constant magnetic field in the direction 61 (i.e., with a constant 
control function S{-) = 6), physical experiments show a translation/rotation of the 
above wall along the nanowire. The corresponding mathematical solution of 
associated with the constant control S, is 




(5) 



where 



/ 1 



\ 



Rg 



COS 6 — sin 9 



\ sin 6* cos 6* / 
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is the rotation of angle 9 around the axis Mei. Furthermore, Equation is invari- 
ant with respect to: 

• translations x t—^ x — a, along the nanowire; 

• rotations Re around the axis ei. 

Hence, denoting A — {9, a), one has a two parameter family of symmetries given 
by M\u(t,x) — Rgu{t,x — a). Therefore, we have a three-parameters family of 
particular solutions of © defined by 

(6) u^^^''^{t,x) = Mf,u\t,x) = Rst+eMo{x + St - a) 
and called travelling wall profiles. 

Theorem 1. There exist eo > and So > such that, for all 61,62 G IR satisfying 
\6i\ < 6q, i = 1,2, for all ai,a2 G M, for every e £ (0,eo), there exist T > and a 
control function 6{-) € L°°{1R,1R) such that, for every solution u of ^ associated 
with the control 6{-) and satisfying 

(7) 301 I ||m(0,-)-w'^'''^'"^(0,-)IIh^ <£, 
there exists a real number 92 such that 

(8) \\u{T,-)-u'^^'^''^^{T,-)\\h2 <e. 

Moreover, there exists real numbers ^2 o,nd a2, with 1^2 ^ ^zl + ^ cr2| < £, such 
that 

(9) \\u{t,-)-u'^^'>'^{t,-)U^ 0. 

t — *4-oo 

The control law 6{-) realizing the conclusion of the theorem is actually given by 
the piecewise constant function 



(10) 6{t) = 



62 ~ if < t < T, 

62 if i > T. 



It combines the advantages of being very simple to implement, and of sharing ro- 
bustness properties in norm, as claimed in the theorem. 

The time T of the theorem is arbitrary, but must be large enough so that 



(72 - CTi 



<6o. 



T 

This theorem shows that the family of travelling wall profiles ^ is approximately 
controllable in norm, locally in 6 and globally in a, in time sufhciently large. The 
controllability property with respect to 9 is not clear. Intuitively the system should 
not be controllable in 9, however this question is not very relevant from the physical 
point of view, since it is the position of the wall which is physically interesting. In 
particular, our result asserts that it is possible to pass approximately (up to the 
variable 9 from a wall profile u^-^-'^ to any other by means of a scalar control of 
the form H10|l . This approximate controllability result may have applications for 
magnetic recording. Note that, on the one part, an exact controllability result does 
not seem to be reachable, due to the physical properties of the system, and on 
the other part, this approximate controllability property is sufficient for practical 
interest. 

Up to now, only the one dimensional case, that is, a ferromagnetic nanowire, 
has been considered for control applications. What happens in the two dimensional 
case is an open question. 
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2. Proof of Theorem 

We follow the same lines as in [3 |H|, and first express the Landau-Lifschitz 
equation in convenient coordinates. This permits to establish stability properties, 
and then to derive the result. 

2.1. Expression of the system in adapted coordinates. The control function 
(|10|l considered here being piecewise constant, it suffices to consider Equation Q 
on each subinterval. Hence, we assume hereafter that the control function S(-) is 
constant, equal to S. Let it be a solution of ||2Jl. Set v{t,x) — R^st{u{t,x — 5t)). It 
is not difficult to check that 

(11) = -u A h{v) - w A (w A h{v)) - 5{vx + viv - ei). 

Consider the mobile frame {AIo{x), Mi{x), M2), where Mi(-) and M2 are defined by 

1 



Mi{x) = I I and M2 = I 1 

-tha;> 



In what follows, we will prove that v is close to AIq. This allows to decompose 

V 



: IR^ X M — > S*^ C IR^ in the mobile frame as 



vit,x) = y/l - ri{t,x)^ - r2{t,x)^Mo{x) + ri{t,x)Mi{x) + r2{t,x)M2. 

Easy but lengthy computations show that v is solution of (|ll|l if and only if r = 
satisfies 

(12) n = Ar + Rs{x,r,rx,rxx), 

where 
(13) 

Rs{x,r,rx,rxx) = S (^q ^) ^ G{r)rxx + Hi{x,r)rx + H2{r){rx,rx) + Ps{x,r) 



and 



• A = (^^j^ ^ j with L^dxx + {1- 2th 2a;)Id; 

• dx + t^xld] 

• G{r) is the matrix defined by 

/ rir2 rl 



G{t) = 



rir2 



+ v/T^M^-i\ 



• Hi(x,r) is the matrix defined by 



^^'2^1 - |f p - rirl -r2+r2r\ 

Hi{x,r) = 



v/1- |r|2ch: 



r2-rl ^1 - |rPr2 + rir^y 



H2{r) is the quadratic form on _ZR defined by 
iJ2(r)(X,X) = pW2 ^ 
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Ps{x,r) = , with 

Vim) 

Pl[.,r) =2..(v/r^- 1)^ - 2.i..|^ - 2.,|.P-A- 



2rl^r^j^ + + r2(l - v/T^f ) + nrl 
5 (^(VT~FF- 1 + rl) + (v/T~M^- l)rith: 



and 



1 ^nshx 



Pl{x,r) = - 2ri(v/l-|r|2-l)— — + 2r2— — - 2r2|ri2 



ch^a; ch^x ch^a 



It is not difficult to prove that there exists a constant C > such that, if 
Ikll|j2 = |rp < i and \d\ < 1, then, for every x & M, for aU p, g € 1?^, 

||i?5(x,r,p,g)||jR2 

^^^^ < c{\5\\\p\\h^ + Ikll^^lklk^ + IklUHblU^ + WAUAMm^ + IIHIrO- 

This a priori estimate shows that Rs(x,r,rx,r,j;x) is a remainder term in Equation 
l)12|l . The rest of the proof reUes on a spectral analysis of the linear operator A, so 
as to establish stability properties for Equation H12() . 

First of all, notice that i is a selfadjoint operator on L^{1R.), of domain H^{1R), 
and that L = with ^ = 9x + thxid (one has £* — ~dx + tha;Id). It follows 
that L is nonpositive, and that kerL — ker£ is the one dimensional subspace of 
L'^{]R) generated by In particular, the operator L, restricted to the subspace 
E = (kerL)-'-, is negative. 

Remark 1. It is obvious that, on the subspace E: 

• the norms \\{-L)^^'^ f\\L^{R) and W/Wh^r) are equivalent; 

• the norms ||L/||i2(jp') and ll/lli/sfjfj-) are equivalent; 

• the norms \\{—L)^^^ f\\L^{iR) and ||/||i/3(jfj) are equivalent. 

Writing A — JL, with 

J ' 



I 1^ 

it is clear that the kernel of A is kerA — kerL x kerL; it is the two dimensional 
space of L'^{IB?) generated by 

ai(x) = and 02(2;) = 

Moreover, combining the facts that L\(^^^j.^± is negative and that Spec J = {1 + 
i, 1 — i}, it follows that the operator A, restricted to the subspace £ = (ker A)^, is 
negative. 

These facts suggest to decompose solutions r of \12\ as the sum of an element of 
ker A and of an element of 8. 
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To this aim, recall that, since Equation (|lll) is invariant with respect to trans- 
lations in X and rotations around the axis ei, for every A = {9, a) e M^, Ma{x) = 
ReMo{x — a) is solution of ((TT|l . Define 



Ra{x) 



({Mk{x),Mi{x)) 
\ (Ma{x),M2) 



the coordinates of M\{x) in the mobile frame (Mi{x), M2(x)). 
We claim that the mapping 

^ :1R^ x£ — > H^{R) 

(A, W) I — > r{x) = Ra{x) + W{x) 

is a diffeomorphism from a neighborhood U of zero in IR^ x £ into a neighborhood 
V of zero in H^{]R). Indeed, if r = R\ + W with W G £, then, by definition, 

(15) (r, ai)i2 = (i?A,ai)L2 and (r, a2)L2 = (i?A, a2)L2 • 

Conversely, if A G J?^ satisfies ((UH)), then W = r - Ra e £. The mapping 
h : IR^ — > M^, defined by h{A) = ((i?A, 01)1.2, (i?A, 02)1.2) is smooth and satisfies 
dh{0) = —2 Id, thus is a local diffeomorphism at (0,0). It follows easily that 5* is a 
local diffeomorphism at zero. 

Therefore, every solution r of H12|) . as long as it stays^ in the neighborhood V, 
can be written as 

(16) r{t,-)=RAit){-) + W{t,-), 

where W{t, ■) € £ and A{t) G 1?^ for every t > 0, and (A(i), W{t, •)) G U. In these 
new coordinates^. Equation (|12l) leads to (see j7j for the details of computations) 

Wt{t,x) ^ AW{t,x)+TZ{6,A{t),x,W{t,x),W,{t,x),W,^{t,x)), 
A'{t)^M{A{t),W{t,-),W^{t,-)), 

where U : R x JR"^ x M x {H'^{IR)f x {H^{R)f x {L'^{R)f — > £ and X : 

R^ X (H^R)) X {L^{R)) — > R^ are nonlinear mappings, for which there exist 
constants > and 77 > such that 



(18) 



\\Tl{5,A,;W,W^,W^^)\\^^H^^m? 

<k{\\A\U- + \5\ + \\W\\ 



(19) \M{A,W,W^)\<k(\\A\\jh2 + \\W\\ 

for every W G £, every 6 £ R, and every A G R^ satisfying ||A||jfj2 < 77. 

Remark 2. Using the fact that L is selfadjoint, it is obvious to prove that AW G £, 
for every W € £; hence, (|17|) makes sense. 



^This a priori estimate will be a consequence of the stability property derived next. 

^This decomposition is actually quite standard and has been used e.g. in |18l to establish 
stability properties of static solutions of semilinear parabolic equations, and in |21 1181 to prove 
stability of travelling waves. 
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2.2. Stability properties, and proof of Theorem ^ We are now in position 

fWi\ 

to establisli stability properties for system (|17|1 . Denoting W ~ [ ), define on 



{H'^{E)f X tlie function 



(20) V{W) = - 



L 
L 



W 



W2 



\\\LW,\\l.^„^^ + \\\LW2\\h(,ay 



Remark 3. It follows from Remark^that, on the subspace E — (kerA)^, yVjW) 
is a norm, which is equivalent to the norm II W^||f//2CjR2 



Consider a solution (VF, A) of |(T7|), such that VK(0, •) 
Since L is selfadjoint, one has 

(21) 
d 



W^o(-) and A(0) = Aq. 



dt 



V{W{t,-)) 



AW, 



''W2 

'(-i)i/2 
{-L) 

Concerning the first term of the right-hand side of (|21|l . one computes 



AW, 



L'^Wi 



' \L''W2 



= -||(-i)^/'W-l||(L2(«))2 - ||(-L)^/2W^2||(L2(K))2 



and, using Remark ^ there exists a constant Ci > such that 



(22) 



, L'^Wi 



<-Cl||iy||^^3(«))2 



Concerning the second term of the right-hand side of (|21() . one deduces from the 
Cauchy-Schwarz inequality, from Remark Q and from the estimate (II8II . that 



(23) 



(-i)V2 







(-L)V2j^(^'A'-'W^'W^^'^^^)'(^;_43/2^, 



(-Lf/^Wi 



< \\n{S,A,;W,w^,w^,)\\ 

(m(R)f l|W^II(if3(iR))2 

< /^(||A||j„2 + |<5| + ||M^||(^2(«)) 

Hence, from (|5lJ, (EH, and (|23l, one gets 



^ll(H3(iR))2 



^V(W^)< + (||A||R2-f|(5| 



IW^I 



\w\ 



If the a priori estimate 



|A(t)|U2-|-|<5| + ||Ty(i,.)||(ff2(K))2 



< 



2K 



holds, then 



j^nwit,-)) < -^||PF(t,.)||^^3(«))2 < -yII^(^'-)II(h2(k))2 < -C2V{wit,-)) 

(using Remark|21for the existence of a constant C2 > 0). It follows that there exist 
constants C3 > and C4 > such that, if \6\ + ||VF(0, •)ll(_f/2(K))2 < |^ are small 
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enough, and if the a priori estimate 

(24) max||A(.)||«.<§ 
holds, then 

(25) \\W{s,-)\\^HHM)r < C3e-^^1|VF(0,-)||(^2(K))2, 

for every s £ [0, T], and moreover, one deduces from l(T7| . ifTi^ . and that, if 
the a priori estimate H24|l holds, then 

||A(t)||^. < ||A(0)||«. + ^||l^(0,.)||(^.(«)).^*e-^^^rfs 

(26) +KCi\\WiO,-)\\lH.^„,)). j\-'^-^ds 

< I|A(o)||r. + ^\\wio, -n^HHmf + ^^11^(0, •)ll(^.(«,))- 

From all previous a priori estimates, we conclude that, if \S\ + ||A(0)||jff2 + 
ll^(Oj ■)ll(_f/2(ji;))2 is small enough, then ||A(f)||R2 remains small, for every t > 0, 
and ||VF(t, •)ll(j?2(jR))2 is exponentially decreasing to 0. 

The first part of the theorem, on the interval [0, T], easily follows from the above 
considerations. For the second part, observe that, from l(T7|l . l(T^ . and (|^ . one 
deduces that ||A'(i)|| jj2 is integrable on [0, +oo), and hence, A{t) has a limit in IR^, 
denoted Aoo = (^oo, foo), as t tends to +oo. The theorem follows with 6*2 =^2+^00 

and a2 = (T2 + (Too ■ 
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